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Abstract. The Goldberger-Miyazawa-Oehme sum rule is used to extract the pion-nucleon coupling constant 
from experimental 7rN information. Chiral perturbation theory is exploited in relating the pionic hydrogen 
s-wave level shift and width results to the appropriate scattering lengths. The deduced value for the 
coupling is f 2 — 0.075 ± 0.002, where the largest source of uncertainty is the determination of the s-wave 
ir~p scattering length from the atomic level shift measurement. 

PACS. 13.75.Gx Pion-baryon interactions - 14.20.Dh Protons and neutrons 



1 Introduction 

The Goldberger-Miyazawa-Oehme [1] sum rule relates the 
coupling constant of the pion-nucleon interaction, / 2 , to 
the difference of the ir~p and n + p s-wave scattering lengths 
and to the weighted integral of the difference of the ir~p 
and tt + p total cross sections. By evaluating the forward 
dispersion relations at the physical threshold, u — fi, we 
obtain [2] 

(1 + »/m){a n - p - a^+p)/^ = 4/ 2 /(> 2 - uj%) + 2 J- ,(1) 

where fi is the (charged) pion mass, m is the proton mass, 
lob = —[i 2 /2m and 

poo 

J- = 1/(4tt 2 ) / (a™-aX)/udk. (2) 
Jo 

The s-wave ir~p and ir + p scattering lengths are denoted 
by a^-p and a^+ p respectively and lu = y fJ- 2 + k 2 , where 
k is the pion laboratory momentum. 

It is well known that eq. (1) is not a precise means to 
determine the value of the pion-nucleon coupling f 2 , but 
it is the natural first step in the process of extracting the 
value from the pion-nucleon data. Furthermore, the sum 
rule makes the connection between the uncertainties in 
the data and in the coupling constant more transparent. 

It is the aim of the present paper to examine the ingre- 
dients of eq. (1), two scattering lengths and integral J~, 
and determine the pion-nucleon coupling strength f 2 . The 
implications of the pionic hydrogen measurements to the 
sum rule are discussed in sect. 2, the s-wave 7r + p scatter- 
ing length is addressed in sect. 3 and the integral J~ in 
sect. 4. Section 5 combines the information to extract a 
value for f 2 and in sect. 6 the conclusions are drawn. 



2 Pionic hydrogen 

The accurate measurements of the properties of pionic hy- 
drogen at PSI provide a source of information on the pion- 
nucleon amplitude very close to the physical threshold. 
The strong interaction level shift ei s , determined through 
the difference of the electromagnetic and measured tran- 
sition energy difference 

has been measured to a precision of about 0.2 % [3] 

e u = -7.120 ± 0.008 ± 0.009 eV. (4) 

The first error is due to statistics and the second due to 
systematics. There is a considerable improvement in ac- 
curacy compared with the earlier result [1] 

e u = -7.108 ± 0.013 ± 0.034 eV, (5) 

mainly because the molecular effects are better under con- 
trol. However, the improvement is not reflected in the ex- 
traction of the n~p s-wave scattering length. The connec- 
tion between the level shift and the scattering lengths is 
provided by the formula [5] 

e ls = -2a 3 n 2 c (a+ + +a Q+ )(l+6 e ), (6) 

where /x c is the reduced mass of the ir~p system, a ~ 
1/137.036 is the fine structure constant and aj + , are 
the isoscalar and isovector s-wave 7rN scattering lengths 
respectively. The quantity S e evaluated next-to-leading or- 
der in isospin breaking and in the low-energy expansion 
has the value [S] 

8 e = (-7.2 ±2.9) x 10~ 2 . (7) 
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The uncertainty of the ir~p scattering length determined 
with eq. (6) is then dominated by the uncertainty in 5 e , 
which in turn is dominated by the largely unknown low- 
energy constant fi [5] . Potential models produce typically 
smaller correction factors S e with considerably smaller un- 
certainties [6l[7]- However, here we rely on the chiral per- 
turbation theory result [5] for S e and consequently the 
measured level shift of eq. (4) gives 



0.0933 ± 0.0029 



(8) 



if errors in eq. (4) are added linearly. Here the identifica- 



tion a 



TT~p 



*0+ 



has been made. 



The Deser-type formula [8] relating the isovector s- 



wave scattering length a 0+ to the level width is [6] 
r ls = 8aVc<7o(l + l/P)[a 0+ (l + S r )} 2 , 



(9) 



where go is the centre-of-mass momentum of the 7r° in 
the charge exchange reaction and P — 1.546 ± 0.009 is 
the Panofsky ratio [9], i.e. the ratio of the cross sections 
a(ir~p — * Tr n)/(r(ir~p — > 771) at the threshold. The cor- 
rection factor Sp has been evaluated in leading order in 
chiral perturbation theory with the result |10j 



S r = (0.6 ±0.2) x 10" 



(10) 



A potential model calculation [6] would give a different 
sign, but this correction factor is rather small in any case. 
The factor Sr depends on the low-energy constant f2 and, 
therefore, the uncertainty is smaller than for the factor S e , 
which depends on the constant f\. 

The s-wave isovector scattering length can be solved 
from eq. (9) 



■'0+ 



0.08933 



1 + 5 r V 0.868 eV 



1/ M . 



(11) 



The measured result for the width [1] _Ti s = 0.868 ± 
0.040 ± 0.038 eV then gives for the scattering length 
a~ + = 0.0888 ± 0.0040 l//x. The relatively large error bar 
is largely due to the uncertainty in estimating the Doppler 
broadening. The new preliminary result for the width [llj 
with improvement in this respect, i~i s = 0.823 ±0.019 eV, 
would give = 0.0865 ± 0.0010 l/>. 



3 Discrete phase shift analysis for n + p 

For the 7r + p interaction a discrete phase shift analysis has 
been performed in the range k — 0.077 — 0.725 GeV/c 
at 77 different momenta [12]. Tromborg corrections [13] 
have been used to extract the hadronic amplitudes from 
the experimental data. Forward scattering constraints [T4] 
have been applied iteratively in the analysis. Additional 
constraints from total inelastic cross sections were used 
for partial wave inelasticities. Special care has been taken 
to incorporate in the analysis the experimental resolu- 
tions and acceptances. In this way only a small number 
of points had to be eliminated due to inappropriate angu- 
lar dependence. The normalizations were allowed to float 
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Fig. 1. The n + p s-wave phase shift 5s3i normalized with 
|flS3i|<7- The asterix denotes the KH80 value. The result of 
Matsinos et al. is shown by the triangle. 



according to experimental information on the systematic 
uncertainty. Furthermore, the behaviour of the zero trajec- 
tories for transversity amplitudes was monitored to keep 
a smooth variation in the shape of observables and the 
extrapolation to the physical threshold was stabilized in 
this manner. Figure 1 shows the phase shift 8szi close to 
the threshold normalized with the product of the absolute 
value of the s-wave scattering length as3i and the centre- 
of-mass momentum q. The result can be parametrized 
with analytical forms in three different q 2 ranges such that 
in the lowest energy range the effective range approxima- 
tion is used. The two dividing momenta have been taken 
as free parameters. It turns out that a two-parameter ef- 
fective range approximation is valid up to q ~ 0.14 GeV/c. 
The result for the scattering length 0531 



0-S31 = a Tr+ 



-0.0764 ±0.0014 1/fi 



(12) 



is very stable. This result completes the discussion of the 
left hand side of eq. (1). The result can be compared with 
the recent figure of Matsinos et al. [15] a^+ p = —0.0751 ± 
0.0039 1/fi based on a fit with a low-energy model and 
with the FA02 solution of the GWU/VPI group [16] a^+ p = 
-0.0911 ± 0.0014 l//i. For the KH80 solution [17] of the 
Karlsruhe group a T+p = -0.1010 ± 0.0040 



4 Integral J 

The integral J~ contains input from a number of sources. 
Results from total cross section measurements are avail- 
able in the range 0.16-640 GeV/c for the n~p scattering 
and 0.16-340 GeV/c for the ir+p scattering [18]. These 
data have been corrected for the electromagnetic effects 
up to k — 0.725 GeV/c using the formalism of Tromborg 
[15] . The corrections have been published only up to 0.655 
GeV/c, but here we employ a smooth extrapolation up to 
0.725 GeV/c. The question of electromagnetic corrections 
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Table 1. Contributions to J (nib) of the different high- 
energy ranges of the laboratory momentum k. 



Input 10-350 GeV/c 350- GeV/c 



Hohler [22] 


0.08786 


0.01787 


Donnachie-Landshoff [23 


0.09968 


0.02514 


Gauron-Nicolescu [24] 


0.10665 


0.02012 


PDG [25] 


0.09587 




Present work 


0.09609 





above the Tromborg range has been discussed in a num- 
ber of articles [T9"ll20) . Here we do not, however, correct the 
data above the limit 0.725 GeV/c, but we estimate the cor- 
responding uncertainty in the value of J~ by making use 
of the prescription in |21j . A smooth curve through the ex- 
perimental points is drawn using the expansion technique 
[22] , which guarantees that forward dispersion relations 
are satisfied. Details can be found elsewhere [14]. 

The integral J~ converges relatively slowly and the 
high-energy behaviour of the cross sections (\.± p plays 
a significant role. Here we consider three different Regge 
type parametrizations for the asymptotic forms, i.e. forms 
which give estimates of the total cross sections beyond 350 
GeV/c: the parametrizations of Hohler [25], Donnachie 
and Landshoff [23] and Gauron and Nicolescu [24]. In the 
range 10 - 350 GeV/c the recent parametrization of the 
Particle Data Group [25] has been employed in addition 
to the Regge forms quoted above and the result from the 
fit with the expansion technique [14] . The results of these 
parametrizations have been summarized in table 1. The 
displayed number of significant figures is there only to help 
in the addition of the contributions of the different mo- 
mentum ranges. The momentum range from the thresh- 
old to 10 GeV/c is split into two at 2.03 GeV/c. For the 
low-energy section from the threshold to 2.03 GeV/c the 
experiment covers only part of the range. From the discus- 
sion in sects. 2 and 3 we can, by invoking isospin invari- 
ance, deduce values for the n ± p total cross sections at the 
physical threshold with the results a™ p = 5.80 ± 0.27 mb 

and cr T 2 t = 1.47 ±0.05 mb. In the range from 0.08 to 0.16 
GeV/c there is information on the angular distributions, 
but no modern total cross section data. For this range 
we have performed a phase shift analysis with elastic n^p 
input and fixed-i constraints. Also, results from the dis- 
crete phase shift analysis for n + p discussed in sect. 3 have 
been used as cross checks. The remaining piece from the 
threshold to 0.08 GeV/c has been covered by a smooth 
interpolation. The contributions to integral J~ can then 
be evaluated with the results displayed in table 2. There, 
in addition, results from the phase shift analyses of the 
Karlsruhe group and the GWU-VPI group are shown. A 
detailed discussion of the contributions of different mo- 
mentum ranges can be found, e.g., in ref. [21] , 

Combining the values from different momentum slices 
yields the value for J~ given in table 3. Also, results from 
some earlier evaluations have been displayed there. To as- 



Table 2. Contributions to J (mb) of the low and interme- 
diate energy ranges of the laboratory momentum k. 



Input 


0-2.03 GeV/c 


2.03-10 GeV/c 


KH80 PU 


-1.27853 


0.10691 


KA84 [26] 


-1.31266 


0.13802 


FA02 [E] 


-1.30213 




Present work 


-1.29757 


0.12046 



Table 3. The values for the integral J (mb). 



Source 


J- (mb) 


Hohler-Kaiser [29] 


-1.06 


Koch [27] 


-1.077 ± 0.047 


Gibbs et al. [28] 


-1.051 ± 0.005 a 


Ericson et al. [21] 


-1.083 ± 0.032 


Present work 


-1.060 ± 0.030 



Statistical error only. 



sign an error bar to our result for J - , we note that the 
statistical uncertainty is quite small, about 0.007 mb. It 
is, however, hard to estimate the systematic uncertainty. 
The systematic normalization uncertainties of the data 
have been taken into account in the fitting procedure [14] . 
Also, there remain discrepancies in the data, see the dis- 
cussion, e.g., in ref. |29| . The Zi-resonance gives a major 
contribution to the J~ -integral and in that range the two 
accurate experiments, Carter et al. [30! and Pedroni et al. 
[31j . differ slightly, but in a systematic manner. The im- 
pact of the difference on the J~ is 0.012 mb. We adopt 
this as the number reflecting experimental systematic ef- 
fects. Furthermore, the forward dispersion relations, which 
have been imposed as contraints, do need a value for the 
pion-nucleon coupling constant as input. Here the pion- 
nuclcon coupling constant has been allowed to vary and 
the effect on the error bar for the J~ integral is 0.001 mb. 
For the asymptotic part, where data do not exist, the av- 
erage of the numbers in table 1 has been taken and the 
corresponding error is chosen such that all the displayed 
values fall within the errors. This gives as the estimate 
for the error 0.004 mb. If in the integrand for J - involv- 
ing crj2* p — cj+p a Coulomb correction of the type [2T] is 

adopted in the range 0.725 - 2.03 GeV/c, the value for J~ 
changes by 0.6 %. Adding up all these errors, we obtain 
the quoted uncertainty 0.030 mb. 



5 Results and discussion 

From eq. (1) the pion-nucleon coupling constant can be 
extracted with the result 

[gC 1 + m^ a *~P ~ a ^ + p)^ ~ J ~M 2 ]- (13) 
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The contribution to f 2 from the term involving the differ- 
ence of the ir~p and ir + p scattering lengths is about 2/3 
and the J~ piece about 1/3. With the scattering lengths 
and the integral J~ evaluated in sects. 2-4, the coupling 
constant becomes 

f = 0.075 ± 0.002. (14) 
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manuscript. One of us (PM) thanks the Magnus Ehrnrooth 
Foundation and the Waldemar von Frenckell Foundation for 
financial support. Support from the Academy of Finland and 
the Russian Academy of Sciences exchange grant is acknowl- 
edged. This work was supported in part by the EU Contract 
MRTN-CT-2006-035482, FLAVIAnet. 



This precision can not compete with the expected accu- 
racy from other methods involving the dispersion relations 
for the B-amplitudes, but we shall use this range of possi- 
ble values as input for the phase shift analysis with fixed-i 
constraints. 

In the discussion above isospin symmetry has been 
used only in fixing the threshold values for the total cross 
sections, a statement valid to 0(p 3 ) in chiral perturbation 
theory. By invoking the isospin invariance we can relate 
a ir-p ~ a -K+p = 2a^~ + , where the isovector s-wave scat- 
tering length is accessible through the atomic width mea- 
surement, eq. (11). The numbers for quoted in sect. 2, 
yield the coupling strength in the range / 2 = 0.076 — 0.077, 
i.e. within the error range quoted above. By making use of 
isospin invariance one can avoid the question of the consis- 
tency of the approach when one is extracting the hadronic 
n~p s-wave scattering length in chiral perturbation the- 
ory and the ir + p scattering length in a discrete phase shift 
analysis exploiting the Tromborg approach for the elec- 
tromagnetic corrections. We are aware of one low-energy 
7rN analysis with a complete treatment of both strong and 
electromagnetic effects to third order in chiral perturba- 
tion theory, the work of Fettes and Meifiner [32]. There, in- 
deed, it is found that nonlinear pion-nucleon-photon cou- 
pling terms generate sizeable effects. However, the com- 
parison in ref. [32j for the low-energy n + p s-wave scatter- 
ing with the Karlsruhe results, which also use Tromborg 
corrections, shows that the effect is relatively small there. 

The potential model [6] corrections 5 e = (—2.1 ±0.5) x 
10~ 2 and 6 r = (-1.3 ±0.5) x 10~ 2 would lead to coupling 
values in the range f 2 = 0.077 — 0.078 in agreement with 
ref. EH- 



6 Conclusions 

The GMO sum rule provides a tool to relate experimen- 
tal information on pion-nucleon scattering from different 
sources. Also, it is relatively straightforward to analyze the 
error propagation. With the input from pionic hydrogen 
and a discrete phase shift analysis for the n + p interac- 
tion we obtain for the pion-nucleon coupling constant the 
value f 2 = 0.075 ±0.002. It turns out, that the largest un- 
certainty in determining a value for the pion-nucleon cou- 
pling constant is currently due to the large uncertainty in 
the correction factor 8 € for extracting the the ir~p s-wave 
scattering length from the ei s hadronic level shift mea- 
surement. Assuming isospin invariance we can make use 
of the isovector s-wave scattering length which will 
be determined within 1 % by the PSI group in the future. 
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